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Abstract 

The recent discoveries of several reliable events of high energy cosmic rays at an 
energy above lO^*^ eV raise questions about their path through the nearby universe. 
The two analyses of, on the one hand, the Haverah Park data set including a limited 
set of further events and, on the other hand, the Akeno data set appear to have an 
inconsistent pattern of arrival directions. Both data sets showed some measure of a 
correlation with the supergalactic plane, the locus of cosmologically nearby galax- 
ies, radio galaxies and clusters of galaxies. In order to be able to interprete such 
findings, we need a reasonable model of the true intergalactic magnetic field and 
then can expect to make further progress on the propagation of energetic charged 
particles. Using recent cosmological simulations of structure formation in the uni- 
verse, we estimate the magnetic fields which correspond to the upper limits in the 
Rotation Measure to distant radio sources. Using the one single direct measurement 
of such a magnetic field, near the Coma cluster, we thus estimate that the magnetic 
field strength in supergalactic sheets and filaments may be in the range of 0.1 to 1 
microgauss. If such strengths are realized inside our Local Supercluster, this opens 
up the possibility to focus charged particles in the direction perpendicular to the 
supergalatic plane, analoguously but in the opposite direction to solar wind modu- 
lation. If focusing exists, it means that for all particles captured into the sheets, the 
dilution with distance d is 1/d instead of l/d?, increasing the cosmic ray flux from 
any source appreciably with respect to three-dimensional expansion. This means in 
effect, that we may see sources to much larger distances than expected sofar. This 
effect is relevant only for energies for which the possible distances are smaller than 
the void scale of the cosmological galaxy distribution, in the range possibly up to 
100 Mpc, but presumably less than this distance. 
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1 Introduction 



The discovery of cosmic rays early this century [|22|, ^ spawned many observations of 
these high energy particles, right up to the recent detection ||20|, ^, [l^ of a significant 
number of particles beyond the energy of 10^° eV; for general reviews of these questions 
see, e.g., [0, |,[T3, 0,0,0. 

While the search into the origin of cosmic rays still awaits the final resolution, there are 
many successful steps that have been taken, from the original Fermi-acceleration process 
ISj p!6[] via the argument that the high energy particles ought to be extragalactic |[Tl| , ^, |12 



to the more recent discoveries already mentioned. 

In this brief discussion we propose to concentrate on the arrival directions of the most 
energetic cosmic rays, and their possible correlation with the matter distribution in the 
nearby universe P, 51, . 



2 Expectations for high energy cosmic rays 

For the source region or in our Galaxy the Larmor radius Vg is given by 

Tg ^ 10 E20 BZI kpc, (1) 

where £"20 is the particle energy in 10^° eV, and -B_5 is the magnetic field in units of 10~^ 
Gauss. 

This means that the Larmor radius is larger than the thickness of the Galactic cosmic 
ray disk (about a kpc; for a review of the Galactic magnetic field see and of the 
order of the size of the source region if radio galaxy hot spots of very high luminosity 
are considered 0. Thus, at such energies, propagation through the Galaxy is nearly in a 
straight line path. 

Other important limitations are obviously losses against photon or magnetic field 
backgrounds, and the time required for acceleration (for a review, see Detailed 
calculations for the propagation have been done by a variety of scientists, e.g., Stanev 
[p7| , and several others. 

For intergalactic space the Larmor radius is conveniently scaled to other units and can 
be written as 

Tg ^ 100 ^20 Mpc, (2) 

where the magnetic field strength is obviously given in units of 10~^ Gauss. This means 
that for the typical upper limits derived from Rotation Measure observations (for a review 
see Kronberg [^), the intergalactic propagation is also in a nearly straight line path. 
Therefore, it is meaningful to ask for the arrival directions on the sky, and whether they 
correlate with any known objects or structures. 
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However, losses in the bath of the microwave background radiation (MBR) hmit the 
distance from which particles can realistically come to less than < 50 to 100 Mpc. In other 
words, integrating over a presumed cosmologically homogeneous source distribution leads 
to a cutoff in the summed contributions from all sources near 5 x 10^^ eV. This is the 
GZK cutoff, named after Greisen, Zatsepin and Kuzmin 



19 



38 



Probable sources at distances < 100 Mpc are distributed in the supergalactic plane 
sheet B 



m 
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which is defined by the Local Supercluster of nearby galaxies (< 
30/z~^Mpc). Therefore arrival directions should cluster toward the supergalactic plane 
from energies, where MBR losses become important. Or, in other words, from near 
5 X 10^^ eV the arrival directions should cluster just as the sources in our neighborhood 
do. However, there should be no clustering below this energy. 



3 Test of prediction 

This prediction was made and explored in various lectures late 1994 and early 1995 g, ||, 
|, |3^ , and has been tested using several data sets. 

First, the Haverah Park Array data set was used, and also combined with a fraction of 
other data available at the time |^6[. In this test the question was asked in the following 
way: Given very limited statistics, we can test whether the arrival directions cluster 
better towards the Galactic plane, or the supergalactic plane. The measure of success 
was the distance of the arrival directions to the reference plane, and the distance for 
homogeneously scrambled data (in order to allow for all selection effects). This test 
showed an effect of a correlation at a level somewhat below 3 sigma. 

Thus, the effect was consistent with the prediction and was visible from 4 x 10^^ eV. 

Second, the Akeno Array data set was used These authors found another measure 
to be better as a test: The distance to the Galactic and supergalactic plane did not yield 
any better result than the analysis of the Haverah Park data, but using pairs of events 
there appeared a tantalizing excess of pairs lying directly on the supergalactic plane 
sheet. Thus, here again the effect was consistent with the prediction, however, using pairs 
of events, from 5 x 10^^ eV. 

Third, as presented at this meeting, the combination of all events and again using 
pairs and triplets of events suggested that the supergalactic plane is the region of origin 
of ultrahigh energy cosmic ray events, for a fraction of maybe 15 - 20 %, or possibly more, 
of all events beyond 5 x 10^^ eV. 



However, as pointed out by Waxman et al. |^, not only is there a seeming inconsis- 
tency between the first two analyses, but there is also a deeper difficulty: When one uses 
the actual distribution of galaxies farther beyond the Local Supercluster as a measure 
of possible source directions and distances, the supergalactic plane is not such a good 
approximation anymore, and so one would not really simply expect a direct straightfor- 
ward correlation. The observed correlation is apparently better with the geometric sheet 
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corresponding naively to the cosmologically nearby galaxy population, rather than with 
the actual large-scale galaxy distribution. 



4 The supergalactic Plane 

Galaxies are distributed in the observed universe in a non-homogeneous pattern, in what 
may loosely be described as a network of filamentary superclusters encompassing voids, 
with typical void scales of (30-100) Mpc |TB[. We live in one of the superclusters, 
namely the Local Supercluster [^. It is a flattened condensation of nearby galaxies 
centered at the Virgo cluster extending to 30h^^ Mpc with a scale height of 5h~^ Mpc. 
It is also connected to several nearby superclusters by filaments of galaxies and clusters. 
This observed structure of supercluster-void networks can now be well simulated in large 
computers and also be interpreted [l^: The large scale structure forms as the result of 
gravitational instability, and then the matter flows into the potential wells, into the sheets, 
filaments, where sheets intersect, and nodes, where filaments intersect. This means that 
there is baryonic accretion flow towards the nonlinear structures in cosmological structure 
formation. The velocity of this accretion flow can be as large as about 1000 km/sec. 



independent of the Hubble constant |30 



As a consequence of the accretion flow, the cosmological magnetic field is expected to 
lie mostly along sheets and filaments. Assuming uniformity, Kronberg et al. had derived 
an upper limit for the true intergalactic magnetic field of B < 10~^ Gauss; 
here Lrev,Mpc is the reversal scale of the magnetic field. Allowing for the correlation of 
magnetic field to cosmic structures we can rederive this limit and it transforms an upper 
limit to B < 10^^^°-^ Gauss along sheets and filaments p9|, |3T|. Interestingly there is 



confirmation of a definitive strength of such a magnetic field in one case, in the plane of 
the Coma/Al367 supercluster p6|. This supercluster is about 90/;,~^Mpc away from the 



edge of the Local Supercluster and they are connected by filaments of galaxies. 



Figure 1. A two-dimensional cut in a simulation of the evolution of the cosmological flow 
in a standard cold dark matter (SCDM) universe with total = 1 and baryonic fij, = 0.06. 
The calculation has been done in a box of (32/i~^Mpc)^ volume, and the plot includes 
a region of (16/i~^Mpc)^ with a thickness of 0.25/i~^Mpc. The first panel shows density 
contours, the second panel shows velocity vectors, and the third panel shows magnetic 
field vectors. In the third panel, the vector length is proportional to the log of magnetic 
field strength. 
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5 Supergalactic Modulation and Confinement 



Suppose we approximate the Local Supercluster as a cosmological sheet bounded by two 
plane parallel accretion shocks. Then in the accretion flow and in the sheet we can again 
write the Larmor radius 

r,^ ^ 0.05 E,g,7 BZI (^) Mpc, (3) 

Vtot 

and notice that the Larmor radius for the highest energy particles is smaller than the 
thickness of the supergalactic sheets for B_q> 0.02 h Eigip±/ptot. This means a rather 
small strength of the magnetic field may be sufficient to contain high energy particles in 
the sheets. 

Then the question arises : Is modulation possible in the accretion flow analogous to 
solar wind modulation of cosmic rays? Of course, if such a modulation were possible, it 
would only pertain to the momentum component perpendicular to the sheet. 

The transport equation for energetic particles with z along ±-direction to the sheet 
can be written as 

d 

= Nv± - K± — N. (4) 

oz 

N ^ No exp[-^^^] (5) 

K± 

Note v± < for z positive, where v± is the accretion velocity perpendicular to the sheet, 
and K,± the transport coefficient perpendicular to the sheet. 

The critical question is then: What numerical value can k± possibly have? The trans- 
port coefficient can be written as (the characteristic length scale) x (the characteristic 
velocity scale). If all the magnetic field is perpendicular to the sheet, then the Larmor 
radius and the speed of particle, c, are the relevant scales, then we have a very large 
K±, and no modulation is effectively possible. But if the magnetic field has large parallel 
components, then convective turbulence is dominant probably. It means that velocity 
and length scales of transport are of the order of the accretion velocity {v±) to the super- 
galactic plane sheet and the scale height of sheet (Hsgp), respectively. Then for particles 
with rg± < Hsgp/2 ~ 2.5 h~^Mpc modulation is possible. Thus, if supergalactic modula- 
tion can exist, then we can have weak confinement along the supergalactic sheet. 

In three-dimensional expansion dilution and a gradual weakening of the particle 
fiux by interaction with the MBR lead at large distances to a total cutoff above 5 x IQ-'^'^ cV. 
If some fraction of the particles originated from the sources in the Local Supercluster are 
confined along the sheet, we have only a two-dimensional 1/d dilution but also interaction 
with the MBR. This effect is clearly relevant only beyond 5 x 10^^ eV, since at lower 
energies we can see sources beyond the Local Superclusters. In other words, for confined 
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particles focusing along the two-dimensional sheet is possibly stronger than original source 
distribution beyond 5 x 10^^ eV. 

We offer this possibility as a potential solution to the conundrum posed by the two 
data analyses about arrival directions in the literature. This is a possible explanation for 
the Haverah Park and Ag results. 

6 Tests 

There are a number of tests that can be done in the next few years, or could be done now: 

• In such a picture we have stronger magnetic fields along sheets/filaments of super- 
clusters, of order 

B < l0"6±°-5 Gauss (6) 

Radio polarization observations of cosmologically distant radio sources (group of 
Phil P. Kronberg) will provide the most stringent check. 

• All present and future events beyond 4 x 10^^ eV from the various arrays Akeno, 
Haverah Park, Fly's Eye, Yakutsk, Volcano Ranch, and in the future Auger should 
be combined to repeat the analysis: At this meeting there was a first report of such 
an attempt, with very interesting results. 

• We need to verify specific source candidates, such as the radio galaxies 3C134, 
NGC315, and M87. We also suggested the accretion shocks around the large scale 
structure ClS db possible candidate 

The future of our attempts to understand the origin of these very high energy particles 
promises to be challenging. 
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